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1.0 


INTRODUCTION 


This document reports changes made to the Lewis Chemical 
Equilibrium Program in order to adapt it for use by the 
Coal Gasification Project. The program was developed by 
the Lewis Research Center for analysis of the chemical 
equilibrium combustion in rocket engines. It can be applied 
directly to the entrained flow coal gasification process. 

The particular problem addressed is the reduction of the 
coal data into a form suitable to the program, since the 
manual process is involved and error prone. A similar prob- 
lem in relating the normal output of the program to parame- 
ters meaningful to the coal gasification process was also 
addressed. 

The revisions included in this writeup act as buffers on 
both sides of the Lewis Chemical Equilibrium Program to 
allow direct processing of raw coal gasification data and 
the production of data meaningful in coal gasification terms. 

This document includes a Program Task Description of the 
capabilities created by the revisions of the original 
program. A list is given for the program elements which 
have been altered. Also included are two sample problems 
and a printout of the original program with and without the 
modifications for processing of coal gasification data. 

The work was performed for the Engineering Systems Branch 
(AH33) of the Computer Services Office under contract 
NAS8-31640. The sponsoring organization was the Coal Gasi- 
fication Task Team (PF15) . 

1.1 References 


This program is self-contained and requires no other docu- 
ments for its use and execution. Included for clarification 
are a memo by David Seymour, "Application of Lewis Chemical 
Equilibrium Computer Program to Coal Gasification" in 
Appendix A (Reference 2) . If further revisions of this pro- 
gram become necessary, consult the "Computer Program for 
Calculation of Complex Chemical Equilibrium Compositions, 
Rocket Performance, Incident and Reflected Shock, and 
Chapman - Jouget Detonation" by Gordon and McBride, 

NASA SP-273 (Appendix B, Reference 1) . 
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2.0 PROGRAM TASK DESCRIPTION 


The Lewis Chemical Equilibrium Program was revised to include 
the following capabilities: 

1. Accept data in the form of coal and other fuel element 
compositions and create reactant data compatible with 
the original program. 

2. Store and edit the coal composition data. 

3. Increase the number of elements within a compound from 
four to five. 

4. Increase the number of reactants in the system from 
15 to 25. 

5. Provide a means of removing certain species from con- 
sideration during coal calculations. 

6. Provide for coal data only, that input temperatures be 
in °F and that printout of pressure, temperature, 
enthalpy and entropy be in English units. 

7. Provide, for coal data only, an additional section of 
output listing. 

8. Allow the use of a new pseudo element, inert carbon. 

9. Provide a store and edit feature for namelist list. 

10. Provide that types of input data, as represented by 
various menu items, be enterable in any order. 

11. Simplify loading the program onto the terminal. 

12. Provide a simple means of creating a copy of the sys- 
tem tape with the data from the previous run. 

The program which provided the base for the revisions was 
"Computer Program for Calculation of Complex Chemical Equilibrium 
Compositions, Rocket Performance, Incident and Reflected Shock, 
and Chapman-Jouget Detonation" by Gordon and McBride, NASA SP-273. 
This program had been previously revised to provide interactive 
input and output via the Marshall Interactive Planning System 
(MIPS) . 

A complete listing of the source code is appended so that the 
changes indicated in the Revisions section may be followed. 
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3.0 REVISIONS OF LEWIS CHEMICAL EQUILIBRIUM PROGRAM 

The current revision is discussed in relation to the changes 
required in various subroutines and how they relate to the 
new capabilities. Also discussed are the new inputs and 
outputs required or generated by the revision. Finally , 
a complete sample case, for coal data, is presented. 

3. 1 Storage Element A (Map overlay source code ) 

A new segment COALCV was added to handle the conversion of 
coal data to reactant data. Note that this segment is the 
longest and may at some future time have to be split into 
more segments. 

3.2 S torage Element BLOCK1 (Miscellaneous data - Block 
Hata form) 


The entry for atomic weight and valence of Deuterium was 
replaced with the atomic weight and valence of the pseudo 
element, inert carbon. The element symbol is IC and the 
weight and valence are those of regular carbon. The inert 
carbon element is used in the simulation of coal gasifier 
systems which have an unreacted carbon recycle stream. 

3 . 3 Storage Element BOOT (B oots trap control section) 

The control statements serve to assign, copy and catalog 
the five files necessary to run the program. The catalog 
is included to insure that additional runs made on the same 
day do not have to reload the tape. 

Note that BOOT is part of the first file on the system 
tape for the program. The following control statements 
are necessary to load the complete system^ 

@ASG ,T PUR , U9S , XXXXX 
3COPY , G PUR.TPF5. 

@ADD TPF$ . BOOT 

The first part of Sample 1 shows the complete loading process. 

File LEWIS contains all the subroutines and map elements. 

File 4 contains the thermodynamic data. Files 8, 10, and 
12 contain dummy reactant data, namelist data and coal 
data, respectively. 

3. 4 Storage Element CHECK (Subroutine CHECK) 

No changes. 
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3. 5 Storage Element COALCV (Subroutine COALCV) 


The subroutine COALCV is a new subroutine designed to do 
the major processing of converting coal composition data 
to reactant data. 

The appended memo "Application of the Lewis Chemical Equili- 
brium Computer Program to Coal Gasification" by David C. 
Seymour (A-l) discusses the detailed calculations to be made 
based on the input data. 

The subroutine assumes the existence of a previous set of 
coal data, stored in file 12. A set of coal reactant data 
immediately follows the coal data on file 12. This data 
is in the same form as the reactant data written by sub- 
routine RREAD for data which is entered in the normal 
fashion. A dummy set of coal data is provided as part of 
the system tape. The user is queried as to whether he 
wants a printout of the previous data set. A response of 
YES produces a listing of the original coal data in the 
form shown in Sample 2. Regardless of whether a printout 
of the past data is requested, the past data is always 
read except for the conditions of first pass and no request 
for past data. 

Table 1 explains the variables used by COAI£V to store the 
coal data and the record groupings of the data. Table 2 
explains the various possible measured or empirical values 
of higher heating value. 

After printing the existing data, the user is prompted as 
to whether he wishes to modify the data, add to the data 
or continue with the same data. For the continue option, 
the subroutine exits. For the modify flag, a further 
prompt is given requesting the data group numbers of the 
data to be changed. See the last line of Sample 1. Each 
data group has associated with it a -flag which indicates 
whether that type of data is to be processed. For the 
modify option, the flags of only those groups requested 
are turned on. For the add option, all flags are turned on. 

Each data group flag is queried. If the flag is on, ap- 
propriate prompts are issued and the resulting information 
as entered by the user is stored in the appropriate varia- 
ble. Two of the data groups, ash and other fuels, have 
additional special processing. For those data groups, 
two paths are possible. If the modify option was chosen, 
then it is possible to modify, delete or add individual 
lines to the existing material, via a series of prompts 
and responses. If the modify option was not chosen, the 
prompts and responses insert a completely new set of data. 
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Note that for ash composition, four of the common compounds are 
prompted for individually. If the weight of any of these is re- 
turned as 0 (or blank), that compound is eliminated. Note also 
that for other fuels, a search is made to determine which fuel 
number corresponds to water so that the proper water to coal 
ratio may be computed. 

When all the required data groups have been entered, then the 
calculations outlined in the included paper are done, except for 
the section "Use and Interpretation of CEC Results”. That sec- 
tion is done in another subroutine. 

The coal data is saved on file 12 followed by the computed re- 
actant data. The coal data is in the form shown in Table 1. 

The reactant data is in the same format as for the original 
program, except that the first record contains the coal to fuel 
ratio (CF) in addition to the number of reactants (NREAC) . The 
subroutine then returns to the calling routine. 
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TABLE 1 


COALCV VARIABLES 


Record 1 


NREAC 

CF 

HHW 

WCR 

PER 

Record 2 

CTEMP 

OTEMP 

Record 3 

PROX ( 3) 

COAL (5,2) 


PER 

IE 

HH ( 8) 


Record 4 

NCOAL 

NASH 

NFUEL 

NOXID 

* 

Record 5 

CCOMP(12,5) 

ACOMF (11,8) 

OFUEL ( 13,8) 

0X10(11,5) 


Number of reactants 

Coal to fuel ratio 

Higher heating value of coal Btu/lb 

Water to coal ratio 

Fraction carbon conversion 


o 

Coal temperature ( F) Q 
Oxidizer Temperature ( F) 


Proximate Analysis 
volatiles, fixed carbon, ash 
Coal composition in form element 
and number in order C, H, N, O, S. 

As above in record 2 

Pointer to higher heating value table 
Measured and empirical higher heating 
value table. 

See table 2 for possible values. 


Number of other coal elements, other 
than C, H, N, 0, S. 

Number of ash compounds 

Number of fuel compounds other than 

coal. 

Number of oxidizer compounds. 


Other coal components in the form, 
weight, percent, state, and compound. 
Ash components in the form, weight 
percent and compound. 

Other fuel components in the form, 
weight, percent, temperature, state 
and compound. 

Oxidizer components in the form, 
weight percent and compound. 


* 

Compound in this section indicates that a chemical compound is 
broken down into its component atoms and number of atoms. That is: 
CaCO^ « Ca 1.0 C 1.0 0 3.0 
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TABLE 2 


Possible Higher Heating Values 
Measured Data 
Dulong 1 

145. 44 (%C) + 620. 28 (%H - l/8(%0)) ♦ 40.5(%S) 

Dulong 2 

146 . 0 (%C) + 620.0(IH-l/8(%0) ) 4- 40.5US) 

CPP2 

144 . 9 ( %C) 610.0 (%H-l/8(%0) ) + 55.5(%S) 

Grummel and Davis 

[ 654. 3 (%H) . T%C . %H - 40 . *sl 

too. - tA.h i + 4J4 - 6 ^| Lt t + tJ 

Coal Conversion System Technical Data Base (CCSTDB) 
146 . 58 (%C) ♦ 568. 78 (IH) + 29.4(«S) 

- 6 .58 (%Ash) - 51.53 (%0 + %N) 

Maximum of 2 through 6 

Minimum of 2 through 6 
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3 . 6 Storage Element CPHS (Subroutine CPHS) 


Changes were made in the COMMON area MISC to accommodate 
25 instead of 15 reactants. The variables changed were 
LLMT, BO, BOP, DATA, NAME (first dimension), ANUM (first 
dimension) PECWT, ENTH, FAZ, RTEMP, DENS, and RMW . In 
addition, the dimension of NAME (second dimension) and 
ANUM (second dimension) were changed to scconanodate 5 
elements within a compound rather than 4. Note that the 
last column of the array is used to store a flag indi- 
cating whether or not the enthalpy of the reactant is to 
be calculated. 

3. 7 Storage Element DETON (Subroutine DETON) 

Sane changes as for CPHS. 

3.8 S torage Element EQLBRM (Subroutine EQLBRM) 

Same changes as for CPHS. 

Also increased temporary storage PROW to 25 from 15. 

Added a temperature condition to the error message which 
states that M 100 iterations did not satisfy the conver- 
gence requirement". 

3.9 Storage Element FROZEN (Subroutine FROZEN ) 

Sine changes as for CPHS. 

3.10 S torage Element GAMEFF (Subroutine GAMEFF) 

No changes . 

3.11 S torage Element GAUSS (Subroutine GAUSS) 
rn changes. 

3.12 S torage Element HCALC (Subroutine HCALC) 

Same changes as for CPHS. 

In addition the temporary storage area NUM was redimen- 
sioned 25 x 6. All references to NAM (N,5) or NAME (N,5) 
were changed to(N,6). All array processing of ANUM was 
increased by 1. 
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3.13 Storage Element LTCPHS (Subroutine LTCPHS) 

No changes . 

3.14 Storage Element MAIN (Subroutine MAIN - system driver) 

Same changes as for CPHS in the data area. Also the addi- 
tion of COMMON area CFUEL and the addition of a data area 
( NC'MTS ) which represent those species which can be elimi- 
nated from considers*. ' when working with coal data. See 

Table 3 for a summary of the species. • The upper limit of 
150 possible species is often exceeded when using coal data. 
Tha optional omission species can be used to remove some of 
the least likely possibilities. 

In the body of the subroutine the value of CF is initialized 
to -1.0. If CF remains as -1.0, then the coal data section 
was never entered since that section- will produce a posi- 
tive value of CF. The CF flag is used in other subroutines 
to initiate action required by the presence of coal data. 

All the entries in namelist, INPT2 , are initialized. 

The menu of possible choices was expanded by addition of 
Item 7, 'start program calculation' and Item 8, 'insert 
or edit coal data.' The wording of Items 1 to 6 were re- 
worked to make them consistent. 

A request for Item 8 produces a call to subroutine COALCV 
and processing returns to the menu list. For menu Item 3, 
the call is made to subroutine RREAD rather than to sub- 
routine REACT. This is to make the data entry independent 
of order. The decision on logical flag CALCH was moved to 
menu Item 7. In the move, references to NAME (N,5) were 
replaced with (N,6) . 

Menu Items 4 and 5 were corrected to properly store the 
species name being placed in the omit and- insert lists. 

Menu J. tern 6, concerning namelist, had an edit feature 
added. A dummy namelist is stored by the system on File 10. 
The user is queried as to whether he wishes to see the 
previous namelist. A response of YES prints the previous 
list and then the user is asked whether he wishes to change 
the namelist. The program either stops for new namelist 
data or returns to the menu list. A response of NO in the 
print namelist question, totally initializes the namelist 
and then stops to input new items. Return is to the menu 
list. 


3-7 


Optional 


AL (S) 

AL(L) 

AL CL2 

ALN(S) 

ALN 

C 

C3 

C4 

C5 

H20 (S) 

03 

SI 

SIC 

SIC2 


TABLE 3 



Omission Species 


SI(S) 

SI(L) 

SICL 

SICL2 

SICL3 

SICL4 

SIH 

SIH4 

SIN 

512 
SI20 
SI2N 

513 


3-8 


Menu Item 7 was added to separate Item 6 from the actual 
start of the program so that Menu Items 1-6 and 8 can be 
input in any order. The new section moves the printout 
of the final namelist from Menu Item 6, the determination 
of logical CALCH from Menu Item 3 and the initialization 
of RHOP, VOL and NT from the original initialization area. 

In addition, this section requests the omission of the 
special coal data omit candidates. 

Some special processing of Menu Item 7 is required by the 
exclusion of coal data. If the fraction of carbon con- 
version (PER) is less than 1, the special inert carbon ele- 
ment must be placed in the insert list. If the logical 
flag ERATIO is false and flag CF is positive then all MIX 
values must be multiplied by CF. MIX supposedly repre- 
sents oxidizer to fuel ratio. For coal data, it is more 
convenient to input oxidizer to fuel ratio. The internally 
calculated coal to fuel ratio (CF) is then used to adjust 
to the proper 0/F ratio. Also all namelist temperatures 
read in°F when coal data is being used, are adjusted to 
°K. 

3.15 Storage Element MATRIX (Subroutine MATRIX) 

Same changes as for CPHS. 

3.16 Storage Element 0UT1 (Subroutine 0UT1 entries 
0 UT2, 0UT3, 0UT 4) 

Same changes as for CPHS. 

In addition, for coal data only, an additional line is 
printed at the end of the 0UT1 printout. The line identi- 
fies oxidizer/coal ratio, water/ coal ratio and fraction 
carbon conversion. 

In the 0UT2 section, for coal data only, pressures are 
printed in lbs per square inches , temperatures in °F , 
enthalpy in Btu/lb, entropy in Btu/(lb *F) , and density in 
lb./cu.ft. The program pauses at the end of 0UT2 to pro- 
duce hardcopy of the output. 

In the 0UT3 section, for coal data only, an extra header 
line is printed. It is exactly the same as the line printed 
by 0UT1 . 

The 0UT4 section is new and corresponds exactly to the sec- 
tion "Use and Interpolation of CEC Results." 

Sample 3 shows the printout created by each section of 
OUTl . 
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3.17 Storage Element PROK (PDP Element PROK) 

No changes . 

3 . 18 Storage Element REACT (Subroutine REACT) 

Same changes as for CPHS. 

Addition of common area CFUEL. Dimension changes of 
temporary variables ANAME, V, LLMTS , SBOP changes to array 
processing to handle 25 instead of 15 reactants. Change 
array processing ANUM and NAME to handle extra element 
within reactant. 

The call to subroutine RREAD was removed and placed in MAIN. 
REACT now creates the final reactant data arrays by adding 
reactants placed on the reactant file (4) by the normal 
reactant input and those coal data reactants which are 
stored on the coal data file (12) . The subroutine uses a 
two pass system. The first pass reads the first record 
from the coal data file if CF is -1. This places necessary 
header data into the system. If any normal reactants exist, 
they are read into the appropriate data areas. When all 
the reactants are read, the coal data reactants are read 
and added to the appropriate data areas. 

3 . 19 S torage Element RKTOUT (Subroutine RKTQUT) 

Addition of COMMON area CFUEL. 

A call to 0UT4 if coal data is present. 

3.20 Storage Element ROCKET (Subroutine ROCKET) 

Same changes as in CPHS. 

Added a prompt for namelist RKTINP. 

3.21 Storage Element RREAD (Subroutine RREAD) 

Temporary arrays NAME, ANUM, PECWT, MOLE, ENTK , FAZ , 

RTEMP , FOX, DENS were adjusted to handle 25 instead of 
15 reactants. NAME and ANUM were also adjusted to handle 
an additional element within the reactant. Array processing 
of these arrays was adjusted accordingly. The COMMON area 
CFUEL was added. 
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The first record of the reactant file was amended to include 
the value of CF. The individual reactant lines were ad- 
justed to isolate the code which indicates that the enthalpy 
should be calculated. Previously, it had been recorded as 
part of the NAME array. 

In addition, the edit section was made two pass. In the 
first pass, the data which is added or edited, is placed on 
the reactant file. Data to be edited also comes from this 
file. On the second pass, if one is requested, it is taken 
from the coal reactant data portion of the coal data file. 
Updates and/or additions are made to this data and the 
changed data is returned to the coal data file behind the 
original coal data. Note that the reactant data from the 
two sources is not added together until subroutine REACT, 
which is part of the startup routine. The setting of flag 
FIRST was removed from subroutine RREAD. 

3.22 Storage Element SAVE (Subroutine SAVE) 

Same changes as for CPHS. Addition of COMMON area CFUEL . 

Creation of an additional header line, for coal data only. 

The line contains oxidizer/coal ratio, water/coal ratio and 
fraction carbon conversion. This header occurs just before 
the summation of atoms in the system.. See Sample 2, labeled 
SAVE. 

3.23 S torage Element SEARCH (Subroutine SEARCH) 

Same changes as for CPHS. Added a specific number for number 
of species in system. Also added a pause at the end of 
the species printout to obtain a hardcopy. 

3.24 Storage Element SHCK (Subroutine SHCK; 


Same changes as for CPHS. 

3.25 S torage Element TAPESAVE (Tape Save Control Section) 

The control statements serve to assign a tape, called 
LEWISTAPE, and copy the five system files into the tape. 
The following statement is necessary to call the tapesave 
system: 


@ADD LEWIS. TAPESAVE 
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3.26 Storage Element THERMP (Subroutine THERMP) 

Same changes as for CPHS. Added COMMON area CFUEL. 

Addition of a pause to create a hardcopy of the material 
prior to the output of THERMP. 

For coal data only, a call to entry OUT4 of subroutine 

OOT1. 

Also, a pause at the end of THERMP to create a hardcopy 
before the next data set is processed. 

3.27 St orage Element VARFMT (Subroutine VARFMT ) 

No changes . 

3.28 T hermodynamic Data (File 4) 

Thermodynamic data for the pseudo element inert carbon was 
added to the data. The data was that for C(S) with the 
exception of the compound name and atom name. 




4.0 SAMPLE PROBLEM 1 

The following pages indicate the procedure to load the 
Lewis Chemical Equilibrium Program and input new coal data 
into the system. 

The data set which is being illustrated is the same as in 
the attached memo (A-l) . 
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5.0 SAMPLE PROBLEM 2 


The following pages which could be a continuation of 
Sample 1, indicate the procedures to be used to edit 
existing coal data and/or existing namelist dat^. 

The data set, after editing will still be the same as in 
the attached memo, but additional namelist options are 
illustrated. 
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APPENDIX A 


APPLICATION OF LEWIS CHEMICAL EQUILIBRIUM 
COMPUTER PROGRAM TO COAL GASIFICATION 

BY 

DAVID SEYMOUR 


This copy of a NASA memo describes the calculations of 
input and output used in the sample problems for conver- 
sion to coal gasification data. 


National Aerpn.mt'Cf and 
Snare Admm'MMt'OP 


NASA 


( 


George C Marshall Space Flight Center 

Marshal Sw t* Flight Center Ai.uvtma 
3581C 


icAr^o* PF15 E^O-llS) 


August 28, 1980 


PF1S/ John P. McCarty 
PF1*>/David C. Seymour 

Application of the Lewis Chemical Equilibrium Computer 
Program to Coal Gasification 

(a) ’’Computer Program for Calculation of Complex Chemical 
Equilibrium Compositions, Rocket Performance, Incident 
and Reflected Shocks, and Chapman - Jouguet Detonations'* 
NASA SP-273, Gordon and McBride 

(b) ’Fortran IV Computer Program for Calculation of Thermo- 
dynamic and Transport Properties of Complex Chemical 
Systems’*, NASA TN-D-7056, Svehla, McBride 

In the early 1960’s Lewis Research Center developed a chemical equilibrium 
combustion (CEC) computer prbpram for rocket engine performance analysis 
(Ref. (a)). Because of the thorough documentation, flexibility, and compre- 
hensive specie library, the program has become an aerospace industry 
standard. 

The CEC program can be applied to entrained flow coal gasification without 
modification. The complete thermodynamic state can be described for either 
adiabatic combustion or specified temperature and pressure. The procedures 
used to apply the program to coal gasification are described here, together 
with a sample case. Because these procedures are somewhat tedious and 
error prone, a subroutine has been added to the CEC program which contains 
these procedures 

The CEC program requires the chemical formula, relative amount and enthalpy 
of each reactant in the system. For a typical entrained flow gasifier, the 
reactants include steam, oxygen (with nitrogen and ar^on impurities) and coal. 
Due to the complex composition of coal, the* chemical formulae for the coal 
components are generally not supplied .md must be determined from an ele- 
mental analysis. 


TO: 

FROM: 

SUBJECT: 

REF: 




A- 1 


Consider the ultimate analysis of the dry Kentucky ^ used by TVA for design 
criteria. 


Coal Composition 


Ash Composition 


Element 

% Wt. 

Compound 

•To wt. 

C 

67. 31 

SiO Z 

45. 94 

H 

4. 757 

A 1203 

17. 31 

N 

1. 529 

F e2°3 

22. 29 

O 

6 . 343 

c a o 

5. 486 

S (Pyritic) 

2. 517 

M g O 

1. 164 

S (Non Pyritic) 

l. 583 

SO 3 

3. 264 

Cl 

. 131 

Na2° 

. 557 

Ash 

15. 83 

k 2 o 

2. 336 


100. 00 

Ti 0 2 

. 742 



Other 

. 911 


. 100.00 


The term ash used here refers to the residue obtained from complete oxidation 
of the coal with excess air during a proximate analysis of the coal. 

The compounds existing in the coal are unknown. In particular, the composition 
of the ash is not the same as the composition of the minerals in the coal and 

thus the ash mass fraction in the; coal indicated above is not the true mineral 

» • 

content of the coal. The sodium and potassium, for example, probably exist 
in the coal as chlorine salts or as part of the other inorganic compounds. Most 
of the iron in the coal exists as iron pyrite, F e S 2 - 

The details of the compounds in the coal and ash are not required for an equili- 
brium ccmbustion calculation. Only the relative amounts of the elements present 
must be known. The ash composition shown above is actually an elemental com- 
position with ficticious oxide compounds assumed and is not a real ash com- 
position. The actual ash composition is dependant on the amount of oxygen avail- 
able during combustion, the combustion temperature, and the rate of cooling of 
the combustion product^. 

The ultimate analysis of coal shown above consists of direct measurement of 
the amount of elemental C, H, N, S, and Cl. The ash weight is then added and 
the oxygen content is taken to be the difference between 100 and this sum. Thus, 
if the orly remaining element in the coal is oxygen, the oxygen content which is 
obtained plus the oxygen in t lie ficticious ash compounds represents the true 
total elemental oxygen in the coal. The ficticious ash compounds ran be used 
directly to describe the relative amounts of the eler nts and corrections to the 
ash content or further assumptions about the coal composition are not required 
The only exception is that the sulfur in the ash should not be included, situ e n 
is already accounted for in the elemental sulfur in the coal. However this 
correction is smaller than the accuracv of t lie data and is not included here. 
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Det e r i nination of CEC Reactant Data F r om Coal Composition 


The elemental mass composition of the coal must be converted into a atom/ 
(mole reactant) form lor input to the CEC. Since the program is limited to 
five elements per reactant, a reacting coal is arbitrarily defined to be com- 
posed of C, H p N p O and S. Other elements and compounds known to exist 
in the coal are treated as individual reactants. 

Since some of the carbon in the coal nir.y not actually be oxidized, a carbon 
conversion factor, f, is defined as the mass fraction of carbon which 
participates in the reaction. The remaining carbon is removed from the 
reacting coal and treated as in inert solid with carbon thermodynamic pro- 
perties. 

As an example of the procedure, the composition of the reacting coal for 
the previously described Kentucky *9 is, with f = . 945 


Coal 

lb Element 

lb Element 

lb Element 

Element 

100 lb ( oal 

100 lb coal 

100 lb reacting coal 

C 

b7. 31 

b3. 61 

79. 18 

H 

4. 757 

4. 757 

5. 921 

N 

1. 529 

1. 529 

1. 903 

O 

6. 343 

6. 343 

7. 896 

S 

4. 10 

4. 10 

5. 103 

Cl 

. 131 

80. 337 

100. oo’ 

Ash 

15.83 




100. 00 


Conve rting 

to lb-atoms / (lb-mole 

reacting coal) gives 



lb Element 

lb-at oms 

lb-atoms 

Element 

lb-atom 

100 lb reacting coal 

100 lb-mole reacting coal 

C 

12. 01 

b. 593 

49. 74 

H 

1. 008 

5. 874 

44. 31 

N 

14. 007 

. 13b 

1. 03 

O 

16. 00 

. 4 

3. 73 

S 

32. 0t>4 

. 1SQ 

1. 20 



13. 25b 

100. 00 


Then the chemical formula feu the reacting coal is 

C 4 q 74 H44 31 N l. 03 °3. 73 S 1 20 

and it's molecular weight per 100 l't-molc is 

(100 ) 100 -= 754. 4 lb / (100 lb-mole) 

13 256 

= 754.4 g(100 g-mole reacting roal) 
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Relative Reactant Weights 


In coal gasification it is convenient to vise as independent variables the mass 
ratios water/coal (or steam/coal) and oxygen/coal. The CEC program, 
however, allows only one independent variable, oxidizer /fuel (O/F). The 
components of the system identified as oxidyzers as opposed to fuels is com- 
pletely arbitrary. The distinction between oxidvzer and fuel is only important 
when running multiple O/F values within a single execution of the program. 

The convention selected here is fuel = coal 4 H 2 0 * (other fuel) , 

compounds 

where coal = reacting coal 4 minerals 4 inert C 4 (other coal ). 

compounds 


Oxidvzer - 0 2 4 (impurities in 0 2 feed), 

where the other fuel and coal compounds and the 02 impurities must be 
included in the CEC thermodynamic data library. 

Then 

coal - J 

fuel 1 4 H 2 O 4 (other fuel compounds ) 
c oa 1 c oa 1 

and t lie relative weight of reacting coal in t lie fuel, required for input 
to t lie CEC is 


reactin g coal = c oal ( reacting c o al) 
fuel fuel coal 

similarly, for t lie other components of the fuel, 

mine ral s = cj^a 1 ( 1 ninerals ) 
fuel fuel coal 

i nc rt C - 1 oal (1-f) (c a rb on) 
fuel fuel < oal 

H 2 0 = coal (H z O) 
fuel file F 1 oal 


and for the oxidvzer 


CEC O/F r cjoaJ (O z ) 

fuel coal 


where 02 /coal includes the impurities. 


As an example, suppose an equilibrium calculation were to be run for 
Kentucky f 9 with f = . 945, H 2 O/ coal = . 10, and 02 /coal - . 9\ t where the 
O 2 feed is 95Tr pure O 2 and S r o Ar by volume. Then the relative weights 
are 

coal = 1 = . 9091 

fuel 1 + . 10 


reacting coal = . 9091 (. 8034) 
fuel 

= . 7304 

minerals = . 9091 (. 1583) 
fue 1 

= . 1439 

inert C = .9091 (1- .945) (.6731) 
fuel 

= . 0337 

H z O = . 9091 (. 10) 

?£r 

= . 0909 

Cl 2 = . 9091 (. 00131) 
fuel 

= . 0012 
1 . 000 


O 2 = . 975 

Oxidyzer * 


Ar = . 025 

oxidyzer 

CEC O/F = . 9091 (. 91) = . 8273 

The relative weights of the ash components are, considering only the 
largest three. 


SjP2 = .1439 ( . 4594 ) = .0773 
fuel .8554 

AI 2 O 3 - . 1439 (. 1731 ) = . 0291 
fuel 7 8554 

Fe203 = . 1439 (. 2229 ) = . 037^ 

fuel . 8554 

. 143 l > 

where the remaining ash components have been arbitrarily spread equally 
between these three. 
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Determination of Coal Enthalpy from HHV 


The CEC program requires an enthalpy lor each reactant when calculating 
adiabatic flame conditions Enthalpies are not required for calculating 
equilibrium conditions at specified temperature and pressure. The enthalpy 
of the reacting coal is the heat of formation ai 26°C plus the enthalpy changes 
from 25°C to the coal temperature. HHV is the heat of combustion at 25°C. 
where the products of combustion are CO 2 , H 2 O ( 1 ), and SO 2 . The HHV of 
coal is either measured in a calorimeter bomb test or jbtained from an 
empirical formula based on the coal ultimate analysis. 

The known HHV CO al can be broken down into the reacting components 
selected for CEC input as 

HHV coal = HHV r c ( reacting co al) 4- HHV C ( inert C) 

coal coal 

The oxides in the ash all have an HHV of zero and do not have to be included 
here. 


The HHV for carbon can be calculated from the heats of formation and the 
definition of HHV. For carbon, 


C 4 0 2 


Z5°C y Cp 2 


HHV = (£hfC 


+ A \>{O z ) - 


(£h f C0 2 ) 


Where ^hfC = heat of formation of carbon, taken here to be amorphous 
carbon. AhfC = 4b80 BTU/lb-mole 

Ah C o 2 = -169293 BTU/lb-mole 

HHV C = (4680 4 169293) J[ = 14486 BTU/lb carbon 

12 . 01 


Then the reacting coal HHV is 

HHV 


r c 


HH\ ( oal * 144 86 ( ine rt C) 

c oal 

r cad i ng co al 
coal 


if. in the previous example, HHV coa j is ui\en as 12141 BTU/lb coal then 
HH v rc - 12141 -J.448M 037) 

. 8034 


r L 2 J 4 ijLlliL 

8034 

=■ 14 446 BTU/lb reacting coal 
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The heat of formation of the reacting coal can be obtained from the 
definition of HHV. The reaction is 

reacting coal + excess C>2 — ► X c C02 + X^ H 2 O (l) + X B SO 2 

where X c# X g are the C, H, and S g-atoms per 100 g-moles of 

reacting coal. 

The N and O in the coal are not included since their heat of formation is 
zero. 

Then HHV rc = A hf rc " ( X c^ h fCC >2 + ^h_A h fH 2 ° + X s ^ h fSC>2^ 

2 

where r - q 4051. 8 cal /g-mole CO2 

^hf H ^o ( 1) = “^8 317. 4 cal / g-mole H 2 0 (l ) 

hfso 2 = “70960 cal /g-mole SC »2 

^hf rr = HHV rc - X c (94051.8) - X h (68317.4) - X s (70960) 

where the units ofAhf r( , required by the CEC are caMlOO g-mole reacting coal). 

# * • 

In the previous example, 

HHV rc = (14446 BTU ) ( cal/g ) (754.4 g ) 

lb reacting coal 1.8 BTU/lb 100 g-mole 

= 6. 0540 x 10° cal/ (100 g-mole) 


and 


Ahr = 6. 0540 x 10 6 - (4b. 74) (94051. 8) - (44. 31) (68317. 4) 
rC 2 

- (1. 20) (70960) 


= 6. 0540 x 10 6 - 6. 2769 x 10 6 


- -222890 cal/(100 g-mole reacting coal) 


If the coal temperature were* 2^0 (77 °F), then the coal enthalpy required 
by the CEC would be Ahf r ^. For temperatures up to about 200°F, the 
reacting coal enthalpy can be estimated from the cmperical equation. 


h - Ah, - . 17S (T - 77) 4 . 002 9V. v (T - 77) 4. 0002S (T - 60) - . 072 

rC 1 r r 


r c 


where T is the coal temperature in °F and W v 1S the percent of volitiles in the 
dry, ash free coal given by proximate analysis. 
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If in the previous example, the proximate analysis were 

Volitiles 37.54 37.54 44.60 

Fixed Carbon 46. 63 46. 63 55. 40 

Ash 15. 83 84.17 100.00 

100. 00 

then W y = 44. 60 and if the coal temperature were 100°F 

hrc -Ahf = . 175 (23) + . 0029 (44. 6) (23) = . 00025 (40) 2 - . 072 
* rc 

= 7. 33 BTU 

lb reacting coal 

= (7.33) 754.4 = 3072 cal/(100 g-mole reacting coal 
1. 8 

and the enthalpy used in the CEC, h rc , would be 

h rc = -222890 + 3072 = -21^820 cal/(l00 g-mole) 

The enthalpv of the remaining reactants will automatically be determined by 
the CEC program from a specified reactant temperature. 

V s e and Interpretation of CEC Results 

The CEC output consists of the thermodynamic state variables and the 
specie mole fractions. For coal gasification several other parameters are 
of interest. 

First, the molecular weight output by the program is an effective value meant 
for use in t lie ideal gas equation of state. This effective value is the mass 
of the system per mole of gas in the system and excludes the moles of solids 
and liquids. To convert to a true ratio of mass /mole for the system, 

M w = (CFC M w ) (1 - !> Xj) 

solid 

+ 

liquids 

where = CEC mole fraction. 

Historically, fuel gas properties are usuullx quoted on a standard cubic 
foot basis, with various conventions taken for standard conditions The 
standard used here is a saturated gas at <»0°F and 30 in of mercury. With 
this standard 385.2 SCF will contain l lb-mole of drv fuel gas. If the standard 
were dry gas at o0°F and 30 in. of mercury, 37 u . 5 SCF would cot. tain 1 lb- 
mole of dry fuel gas. 
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The moles of dry product pas per total moles of reactants is taken to be 

X D p = X C 0 4 X H 2 + X CH4 4 X C02 4 X N2 4 x Ar 

The volume of clean, dry pas produced per lb of coal is then 

Ynp r 385. 2 X pp Reactants , SC F/ (lb coal) 

M w coal 

where the mass fraction ( Reactants ) can be obtained as ( fuel + O 2 ). 

coal coal coal 

The higher heating value of the gas, in BTU/SCF, is obtained from the 
HHV of the constituents as 

HHV r 319. 24 X }{2 + 316. 06 X co + 994. 4 X CH4 

X DP 

The lower heating value, which takes water vapor rather than liquid 
water as the reference state, is given by 

LHV = 269. 71 XH2 + 316. 06 XqO + 895. 3 X CH4 

X DP 

» • 

Cold gas efficiency, which is the fraction of the coal chemical energy 
converted to a gas chemical energy, is 

E = HHV V nP 
HHV coal 

where HHV coa i is the higher heating value of dry coal and does not 
include any of the previously discussed corrections for carbon conversion. 

The equilibrium thermochemical data provided by the program represent 
an idealized state and must be applied to actual gasifier operation with 
caution. The major non -equilibr ium effect, unreacted iarbon, had been 
accounted for in this application of the CEC program. However, as the 
gas/particle mixture is cooled within the gasifier, the chemical reaction 
rates will decrease and chemical equilibrium will no longer be maintained . 

A freeze point emperature can be defined as that temperature at which the 
equilibrium composition approximates the final non-equilibrium temperature. 
This freeze point temperature must be determined by a finite rate analysis or 
experimental data. When the freeze point temperature is known, the CPX 
program can be applied with a specified temperature and pressure to describe 
the final gas composition. 
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It transport property data is required, the CKC program described in Ret 
can be used in lieu of Ref. (A). 


David C. Seymour 
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SOD ROE LlSTlNd OF PROGRAM ELEMENTS 
WITH MODIFICATIONS 


Ihis is a printout of the Lewis Chemical Equilibrium Program 
with the modifications necessary for coal gasification data. 
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APPENDIX C 

SOURCE LISTING OF PROGRAM ELEMENTS 
WITHOUT MODIFICATIONS 


This is a printout of the Lewis Chemical Equilibrium Pro- 
gram without the modifications necessary for coal gasifi- 
cation data. There have been intermediate changes, 
principally to make the program an interactive one. 
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